use "/Users/jg74853/Desktop/ZB - Confidence Book/Data/Final Wave 1 Survey/Wave 1 Weighted (Pooled).dta"

***Defining Variables***
replace q8 =. if q8 > 4

generate q8a = .
replace q8a = 0
replace q8a = 1 if q8 < 3

label define CONFIDENCE 0 "Not Confident" 1 "Confidence"
label values q8a CONFIDENCE

generate pid3 =.
replace pid3 = 1 if partyid7 < 4
replace pid3 = 2 if partyid7 == 4
replace pid3 = 3 if partyid7 > 4

label define PARTY 1 "Democrats" 2 "Independents" 3 "Republicans"
label values pid3 PARTY

generate dem = 0 if pid3 < 4
replace dem = 1 if pid3 == 1
generate rep = 0 if pid3 < 4
replace rep = 1 if pid3 == 3
generate ind = 0 if pid3 < 4
replace ind = 1 if pid3 == 2

generate ideo3 =.
replace ideo3 = 1 if ideo < 4
replace ideo3 = 2 if ideo == 4
replace ideo3 = 3 if ideo > 4 & ideo < 8

generate ideo3a = ideo3
replace ideo3a = 4 if ideo == 8

gen lib = 0
replace lib = 1 if ideo3 == 1
gen mod = 0
replace mod = 1 if ideo3 == 2
gen con = 0
replace con = 1 if ideo3 == 3

generate woman = 0
replace woman = 1 if gender == 2
generate male = 0
replace male = 1 if gender == 1

gen vet = veteran if veteran < 98
replace vet = 0 if veteran == 2
replace vet = 0 if veteran2 == 1
replace vet =. if veteran2 == 98

gen activeduty = 0 if veteran < 98
replace activeduty =. if veteran2 == 98
replace activeduty = 1 if veteran2 == 1

generate milvet= 0
replace milvet = 1 if vet == 1
replace milvet = 2 if activeduty == 1

generate post911vet = 0
replace post911vet = 1 if milvet > 0 & vet_date1 < 2001
replace post911vet = 2 if milvet > 0 & vet_date1 > 2000

generate social = 0 if q45 < 98
replace social = 1 if q45 == 1

generate family = 0 if q46 < 98
replace family = 1 if q46 == 1
replace family = . if q46b_dk == 1

generate religion = 0 if relig < 98
replace religion = 1 if relig == 1
replace religion = 2 if relig == 2
replace religion = 3 if relig > 8 & relig < 12
replace religion = 1 if relig == 12

label define RELIGION2 0 "Other" 1 "Christian" 2 "Catholic" 3 "None"
label values religion RELIGION2

generate christian = 0
replace christian = 1 if religion == 1
generate catholic = 0
replace catholic = 1 if religion == 2
generate norelig = 0
replace norelig = 1 if religion == 3

generate white = 0
replace white = 1 if racethnicity == 1
generate black = 0
replace black = 1 if racethnicity == 2
replace black = 1 if racethnicity == 3
generate hispanic = 0
replace hispanic = 1 if racethnicity == 4
generate asian = 0
replace asian = 1 if racethnicity == 6

gen race2 = racethnicity
replace race2 =. if racethnicity == 3
replace race2 =. if racethnicity == 5

label define RACE 1 "White" 2 "Black" 4 "Hispanic" 6 "Asian"
label values race2 racethnicity RACE

generate generation =.
replace generation = 1 if age > 73
replace generation = 2 if age > 54 & age < 74
replace generation = 3 if age > 38 & age < 55
replace generation = 4 if age > 22 & age < 39
replace generation = 5 if age < 23

generate silent = 0
replace silent = 1 if generation == 1
generate boomer = 0
replace boomer = 1 if generation == 2
generate genx = 0
replace genx = 1 if generation == 3
generate milen = 0
replace milen = 1 if generation == 4
generate genz = 0
replace genz = 1 if generation == 5

generate northeast = 0
replace northeast = 1 if region4 == 1
generate midwest = 0
replace midwest = 1 if region4 == 2
generate south = 0
replace south = 1 if region4 == 3 
generate west = 0
replace west = 1 if region4 == 4 

generate city = 0
replace city = 1 if urban3 == 1
generate suburban = 0
replace suburban = 1 if urban3 == 2
generate rural = 0
replace rural = 1 if urban3 == 3

generate unemployed = 0 if employ < 3
replace unemployed = 1 if employ > 2
replace unemployed = 0 if employ == 5

generate income5 =.
replace income5 = 1 if income < 5
replace income5 = 2 if income > 4 & income < 9
replace income5 = 3 if income > 8 & income < 12
replace income5 = 4 if income > 11 & income < 15
replace income5 = 5 if income > 14

gen married =.
replace married = 1 if marital == 1
replace married = 0 if marital > 1

replace q14=. if q14 > 76
replace age=. if age == 286

generate q1aCOR = 0 if q1a < 77
replace q1aCOR = 1 if q1a == 3
generate q1bCOR = 0 if q1b < 77
replace q1bCOR = 1 if q1b == 4
generate q2aCOR = 0 if q2a < 77
replace q2aCOR = 1 if q2a == 3
generate q2bCOR = 0  if rnd_02 == 1
replace q2bCOR = 1 if q2b_1 == 1 & q2b_2 == 1 & q2b_4 == 1 & q2b_5 == 1 
generate q3COR = 0 if q3 < 98
replace q3COR = 1 if q3 == 3
generate q4COR = 0 if q4 < 98
replace q4COR = 1 if q4 == 1
generate q5COR = 0  if q5 < 98
replace q5COR = 1 if q5 == 2
generate q6COR = 0  if q6 < 98
replace q6COR = 1 if q6 == 3

generate q1aCOR2 = q1aCOR
replace q1aCOR2 = 0 if q1aCOR == .
replace q1aCOR2 = . if q1aCOR == 98

generate q1bCOR2 = q1bCOR
replace q1bCOR2 = 0 if q1bCOR == .
replace q1bCOR2 = . if q1bCOR == 98

generate q2aCOR2 = q2aCOR
replace q2aCOR2 = 0 if q2aCOR == .
replace q2aCOR2 = . if q2aCOR == 98

generate q2bCOR2 = q2bCOR
replace q2bCOR2 = 0 if q2bCOR == .
replace q2bCOR2 = . if q2bCOR == 98

generate q2bCORa = .
replace q2bCORa = q2b_1
replace q2bCORa = 1 if q2b_1 == 0 & q2b_2 == 1
replace q2bCORa = 1 if q2b_1 == 0 & q2b_2 == 0 & q2b_5 == 1

generate q2bCORb = .
replace q2bCORb = 0 if q2b_1 < 98
replace q2bCORb = 1 if q2b_1 == 1 & q2b_2 == 1
replace q2bCORb = 1 if q2b_1 == 1 & q2b_5 == 1
replace q2bCORb = 1 if q2b_2 == 1 & q2b_5 == 1

generate knowledge = (q1aCOR2 + q1bCOR2 + q2aCOR2 + q2bCOR2 + q3COR + q4COR + q5COR + q6COR)
gen k2 = knowledge/6

*****Reversed for Table*****
replace q18 =. if q18 == 98

gen q18a = 0 if q18 < 77
replace q18a = 1 if q18 > 2 & q18 < 77

gen d2 = 0
gen d3 = 0
gen d4 = 0
gen d5 = 0
gen d6 = 0
gen d7 = 0
gen d8 = 0

replace d2 = 1 if dov_assignment_a == 2
replace d3 = 1 if dov_assignment_a == 3
replace d4 = 1 if dov_assignment_a == 4
replace d5 = 1 if dov_assignment_a == 5
replace d6 = 1 if dov_assignment_a == 6
replace d7 = 1 if dov_assignment_a == 7
replace d8 = 1 if dov_assignment_a == 8

label define LOOKS 1 "Mil Looks Like Me" 0 "Doesn't Look Like Me"
label values q18a LOOKS

gen educ5 = educ4
replace educ5 = 5 if educ > 12 & educ < 15

xi I.dem*educ5 I.rep*educ5 I.dem*k2 I.rep*k2, prefix(_I)

gen q16a3 = q16a if q16a < 4

gen q16aR = 0 if q16a3 < 4
replace q16aR = 1 if q16a3 == 2

gen q16aREP = q16aR
replace q16aREP = 0 if q16b == 1
replace q16aREP = 0 if q16b == 3
replace q16aREP = 1 if q16b == 2

replace q44 =. if q44 == 98
replace q45 =. if q45 == 98
replace q46 =. if q46 == 98

gen Q46 = .
replace Q46 = 1 if q46a_1 == 1
replace Q46 = 2 if q46a_2 == 1
replace Q46 = 3 if q46a_3 == 1
replace Q46 = 4 if q46a_4 == 1
replace Q46 = 5 if q46a_5 == 1

gen Q46FAM = 0 if q46 < 98
replace Q46FAM = 1 if q46b_1 == 1
replace Q46FAM = 2 if q46b_2 == 1
replace Q46FAM = 3 if q46b_3 == 1
replace Q46FAM = . if q46b_dk == 77

gen Q46FAM1 = 0 if Q46FAM < 4
replace Q46FAM1 = 1 if Q46 == 1

gen Q46FAMpre = 0 if Q46FAM < 4
replace Q46FAMpre = 1 if Q46FAM == 1
gen Q46FAMpost = 0 if Q46FAM < 4
replace Q46FAMpost = 1 if Q46FAM == 2
gen Q46FAMboth = 0 if Q46FAM < 4
replace Q46FAMboth = 0 if Q46FAM == 2
replace Q46FAMboth = 1 if Q46FAM == 3

gen Q46FAM2pre = 0 if Q46FAM < 4
replace Q46FAM2pre = 1 if Q46FAM == 1
gen Q46FAM2post = 0 if Q46FAM < 4
replace Q46FAM2post = 1 if Q46FAMpost == 1
replace Q46FAM2post = 1 if Q46FAMboth == 1

replace dov_institution = 0 if dov_assignment_a == 1 in 4578/9155
replace dov_institution = . if dov_assignment_a > 1 in 4578/9155
label values dov_institution INSTITUTION

label define INSTITUTION 1 "SCOTUS" 2 "Congress" 3 "Media" 4 "FP Elite" 0 "Military"

gen q15a2 = 0 if q15a < 77
replace q15a2 = 1 if q15a < 3
gen q15b2 = 0 if q15b < 77
replace q15b2 = 1 if q15b < 3
gen q15c2 = 0 if q15c < 77
replace q15c2 = 1 if q15c < 3
gen q15d2 = 0 if q15d < 77
replace q15d2 = 1 if q15d < 3
gen q15e2 = 0 if q15e < 77
replace q15e2 = 1 if q15e < 3
gen q15f2 = 0 if q15f < 77
replace q15f2 = 1 if q15f < 3
gen q15g2 = 0 if q15g < 77
replace q15g2 = 1 if q15g < 3

generate Q13T =.
replace Q13T = 0 if q13 > 2
replace Q13T = 1 if q13 < 3
replace Q13T = . if dov_institution ==.

label define Confidence 0 "Not Confident" 1 "Confident"
label values Q13T Confidence
********************

********************
*****Table 4.1*****
by dov_institution, sort: tab pid3 Q13T [aweight=weight], row

*****Figure 4.5*****
graph bar (mean) q15a2 [pweight = weight], over(dov_institution, label(angle(vertical) labsize(vsmall))) bar(1, fcolor(gs6)) blabel(bar, color(black) size(vsmall) format(%3.2f)) ytitle(% Too Partisan) ytitle(, size(small)) yscale(range(0 1)) ylabel(#6) name(party, replace) graphregion(fcolor(white) lcolor(white) lwidth(thin) lalign(outside) ifcolor(none) ilcolor(none)) plotregion(fcolor(white) lcolor(black) lpattern(solid) lalign(outside) ifcolor(none))

graph bar (mean) q15b2 [pweight = weight], over(dov_institution, label(angle(vertical) labsize(vsmall))) bar(1, fcolor(gs6)) blabel(bar, color(black) size(vsmall) format(%3.2f)) ytitle(% Good at Job) ytitle(, size(small)) yscale(range(0 1)) ylabel(#6) name(good, replace) graphregion(fcolor(white) lcolor(white) lwidth(thin) lalign(outside) ifcolor(none) ilcolor(none)) plotregion(fcolor(white) lcolor(black) lpattern(solid) lalign(outside) ifcolor(none))

graph bar (mean) q15c2 [pweight = weight], over(dov_institution, label(angle(vertical) labsize(vsmall))) bar(1, fcolor(gs6)) blabel(bar, color(black) size(vsmall) format(%3.2f)) ytitle(% Representative of America) ytitle(, size(small)) yscale(range(0 1)) ylabel(#6) name(diverse, replace) graphregion(fcolor(white) lcolor(white) lwidth(thin) lalign(outside) ifcolor(none) ilcolor(none)) plotregion(fcolor(white) lcolor(black) lpattern(solid) lalign(outside) ifcolor(none))

graph bar (mean) q15d2 [pweight = weight], over(dov_institution, label(angle(vertical) labsize(vsmall))) bar(1, fcolor(gs6)) blabel(bar, color(black) size(vsmall) format(%3.2f)) ytitle(% Highest Professional Ethics) ytitle(, size(small)) yscale(range(0 1)) ylabel(#6) name(ethical, replace) graphregion(fcolor(white) lcolor(white) lwidth(thin) lalign(outside) ifcolor(none) ilcolor(none)) plotregion(fcolor(white) lcolor(black) lpattern(solid) lalign(outside) ifcolor(none))

graph bar (mean) q15e2 [pweight = weight], over(dov_institution, label(angle(vertical) labsize(vsmall))) bar(1, fcolor(gs6)) blabel(bar, color(black) size(vsmall) format(%3.2f)) ytitle(% Truthful) ytitle(, size(small)) yscale(range(0 1)) ylabel(#6) name(truthful, replace) graphregion(fcolor(white) lcolor(white) lwidth(thin) lalign(outside) ifcolor(none) ilcolor(none)) plotregion(fcolor(white) lcolor(black) lpattern(solid) lalign(outside) ifcolor(none))

graph bar (mean) q15f2 [pweight = weight], over(dov_institution, label(angle(vertical) labsize(vsmall))) bar(1, fcolor(gs6)) blabel(bar, color(black) size(vsmall) format(%3.2f)) ytitle(% Shares My Values) ytitle(, size(small)) yscale(range(0 1)) ylabel(#6) name(shares_values, replace) graphregion(fcolor(white) lcolor(white) lwidth(thin) lalign(outside) ifcolor(none) ilcolor(none)) plotregion(fcolor(white) lcolor(black) lpattern(solid) lalign(outside) ifcolor(none))

graph bar (mean) q15g2 [pweight = weight], over(dov_institution, label(angle(vertical) labsize(vsmall))) bar(1, fcolor(gs6)) blabel(bar, color(black) size(vsmall) format(%3.2f)) ytitle(% Others Have Confidence) ytitle(, size(small)) yscale(range(0 1)) ylabel(#6) name(respected, replace) graphregion(fcolor(white) lcolor(white) lwidth(thin) lalign(outside) ifcolor(none) ilcolor(none)) plotregion(fcolor(white) lcolor(black) lpattern(solid) lalign(outside) ifcolor(none))

graph combine party good diverse ethical truthful shares_values, name(first, replace)

*****Figure 4.7*****
graph combine respected, name(second, replace)

*****Figure 4.9*****
logit Q13T q15a2 q15b2 q15c2 q15d2 q15e2 q15f2 q15g2 q18a dem rep ideo3 male activeduty vet family social catholic christian norelig white black hispanic asian educ5 boomer genx milen genz midwest south west city rural unemployed income5 married [pweight = weight]
coefplot, drop (_cons) xline(0) xscale(range(-2 2)) xlabel(#5) name (fig49a, replace)

logit Q13T q15a2 q15b2 q15c2 q15d2 q15e2 q15f2 q15g2 dem rep ideo3 male activeduty vet family social catholic christian norelig white black hispanic asian educ5 boomer genx milen genz midwest south west city rural unemployed income5 married [pweight = weight]
coefplot, drop (_cons) xline(0) xscale(range(-2 2)) xlabel(#5) name (fig49A, replace)

clear

*****Run Wave 2*****
